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Abstract 

We consider a situation involving a three-dimensional marginal 
separation, where a (steady) boundary- layer flow is on the verge of 
separating at a point (located along a line of symme t ry/cent re 1 i ne) . 

At this point we include a ‘‘triple-deck’’, thereby permitting a small 
amount of interaction to occur. Unsteadiness is included within this 
interaction region through some external means. It is shown that the 
problem reduces to the solution of a non-linear, unsteady, partial- 
integro system, which is solved numerically by means of time-marching 
together with a pseudo-spectral method spatially. A number of solutions to 
this system are presented which strongly suggest a breakdown of this system 
may occur, at a finite spatial position, at a finite time. The structure 
and details of this breakdown are then described. 


*Work funded by Space Act Agreement C-99066-G. 


In a recent paper, Duck (1989) (hereafter refered to as I), studied 
the effect of allowing a small amount of three-dimensional interaction 
(based on triple-deck theory) at a point where a steady laminar boundary 
layer was on the merge of separating. It was found that this interaction 
zone (i) permitted a small amount of reversed flow to occur within the 
interaction region and (ii) illustrated that non-uniqueness of solution 
was possible. These observations for the three-dimensional case were in 
line with related two-dimensional work published some years previously. 

The original work in this area was that of Stewartson, Smith and Kaups 
(1982), who analysed the (two-dimensional) problem; this was 
studied in further detail by Brown and Stewartson (1983), who showed that 
up to four solutions were possible for a particular choice of parameters. 

This two-dimensional work was later extended to a three-dimensional 
situation, along a line of symmetry by Brown (1985), although, (as 
described in I), a number of assumptions were made to simplify the 
computations, rendering the problem basically two-dimensional. 

The work of Stewartson et al (1982) was extended to the unsteady 
regime by Smith (1982), who showed how the interacting flow may breakdown, 
with the development of a singularity at a point, at a finite time; this 
singularity was also described analytically (in addition to being 
observed numerically). Related studies of problems of this kind were 
then made by Ryzhov and Smith (1984), Smith and Elliott (1985), and extended 
to a rather more non-linear regime by Elliott and Smith (1987). These 
studies all confirmed the possibility of a finite-time breakdown of the 
solution. Further Goldstein, Leib and Cowley (1982) showed how 
two-dimensional Tol lmien-Schl icht ing waves may be generated on marginally 
separated flow. 

A number of other two-dimensional studies indicating the possible 
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occurence of a finite-time singularity of a boundary layer in a triple- 
deck/ inte rac t ing boundary layer include the studies of Smith and Bodonyi 
(1985), Duck (1985a), Duck (1985b), Duck (1987), Tutty and Cowley (1986). 
Recently Brothe r ton-Radc 1 i f fe and Smith (1987) and Smith (1988) have 
presented possible structures of such breakdowns. Further, Duck (1990) 
has presented a number of numerical results indicating the occurence 
of a finite-time breakdown in a three-dimensional non-linear, unsteady 
triple-deck flow. 

In this paper we extend the three-dimensional work of I into the 
unsteady regime. As in I, we concern ourselves with a steady laminar 
boundary layer, which is on the verge of separating at a point, situated 
along a line of symmetry; in particular we assume that both the streamwise 
wall shear and the crossflow derivative of crossflow wail shear (directions 
defined relative to the freestream) vanish simultaneously at the point 
in question, but then both immedially recover. This boundary layer could, 
for example, be set up on a body of revolution, as for example described 
by Cebeci, Khattab and Stewartson (1980) and Cebeci and Su ( 1988). We 
then assume that this flow is perturbed in some (unsteady) fashion. 
Physically, this could perhaps be realised by some small amplitude 
pitching or yawing motion of the body (for example), although for 
this paper the particular details are not important. 

We fully expect (guided by the two-dimensional results of Smith 1982) 
that the unsteady three-dimensional marginal separation problem may 
suffer finite-time breakdowns. The attractive feature of the problem is 
that this breakdown is likely to be partly analysable, and indeed, this 
turns out to be the case; the author is unaware of any previous 
descriptions of a fully three-dimensional unsteady breakdown of an 
interacting boundary- 1 aye r problem. 

The structure of the problem and the numerical techniques employed arc 
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broadly based on I, although the added dimension of time adds significantly 
to the computational complexity of the problem. 

This particular problem is likely of practical importance, being 
linked to the problem of dyanamic stall on an aerofoil or on a body of 
revolution (for example), involving an abrupt change in the 
characteristics of the flow at some finite angle of incidence (for 
example) . 

In the following section we go on to formulate the problem in detail. 
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2. Formula Lion 


Details here are similar in many respects to I, and consequently 
discussion of features similar in this study will be kept to a minimum. 

We consider an incompressible fluid of kinematic viscosity v. 

We take L to denote a typical lengthscale and U„ a typical frecstrcam 
velocity directed along the line of symmetry of the configuration. 
Cartesian coordinates L(x,y,z) are chosen, with origin at the separation 
point; the velocity is then written as U^fu.v.w), and the pressure is 
PoofL 2 p, p„o being the fluid density. The line of symmetry corresponds 
to z = 0, and the body surface (assuming insignificant curvature) is 
taken to lie along y = 0. 

The Reynolds number is defined as 

R = UooL/v, (2.1) 

and this is taken to be large throughout this paper. 

If a (typical) timescale of any unsteadiness is taken to be T, then 
we can define a second non-dimensional parameter, namely 

S 0 = L/TUoo, (2.2) 

and the non-dimensional equations governing the flow may be written as 
V. u = 0, (2.3) 

9u 

So + (u.V)u = - Vp + R'ly2u, (2.4) 

Symmetry of the flow about z = 0 implies u.v.p are even functions 
about this plane, whilst w must be an odd function of z. Throughout, 
on y = 0, we impose the zero velocity condition u = 0. 

The details of the non-interaction zone are precisely the same as in 
I. If we define (standard) boundary- 1 ayer variables Y and V as follows 
y = R-i Y, v = R-i V, (2.5) 

then using the results found in I, as x 2 + z 2 -♦ 0, for Y = 0(1) 
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u = U 0 (Y) + pr U 0 ’(Y) - x[wj (Y)+Vj ’ (Y) ] 

- i x2[V 2 ’(Y) + w 2 (Y) ] 

+ i p 2 r 2 Uo”(Y) + z 2 u 23 (Y) 

- pxr [ V, * ’ (Y) + wfCY)] + (2.6a) 

V = V j (Y) + p ^ U 0 (Y) 

+ x V 2 (Y) - p 2 X 2 x U 0 ’(Y) 

-Hr [ W1 (Y) - Vj ’ (-Y) ] 



+ x 2 ^ [2wj (Y) + V! ’ 

(Y) ] + .... 

(2.6b) 

w 

= z wj (Y) + zx 

w 2 (Y) + 

pzr wj ’ (Y) + . . . , 

(2.6c) 

p 

= P 0 + Pix + i 

P 2 x 2 + 

) P 3 z 2 + . . . 

(2.7) 

Here we 

have wr i t ten 




r 

= (X 2 x 2 +z 2 )i , 



(2.8) 


and we must have that 

U o (0) = U o ’(0) = V j (0) = Vj’(0) 

= wj (0) = wj ’(0) = Vj ’ ’(0) 


» w 2 (0) = V 2 (0) = 0, 

(2.9) 

V 2 ' (0) = P 2 X 2 U o ”(0), 


u 2 3(0) = - i P 2 U o ”(0). 

(2.10) 

Further details are given in I 

; p and X may be regarded 


parameters to the problem. The key outcome of the above equations is the 
existence of a discontinuity of derivatives of the solution as 
x 2 + z 2 -» 0. In order to alleviate this, we introduce a small interactive 
zone, wherein 

X = x/5, Z = z/5, (2.11) 

are the important streamwise and crossflow coordinates respectively, where 

5=R-l/5. (2.12) 

At this stage we must be specific about the magnitude of the 
unsteadiness parameter S 0 . If unsteadiness is to play an important 
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role in the problem (i.e. a time derivative is present in the key 
equation (2.31), derived below), we must have that 

S 0 = 5 1 / 4 S, S = 0(1), (2.13) 

representing a relatively long timescale. A similar scaling was found 
by Smith (1982) to be important. 

For Y = 0(1), the solution in the interaction zone is expected to 
develop as 

u = U 0 (Y) + 5 { - X[wj(Y) + V! * (Y) ] 

+ A(X.Z.t) U 0 ’ (Y) | 

+ 5 2 { - $X 2 [ V 2 ’ (Y) + w 2 (Y) ] 

+ | [A(X ,Z , t ) 2 + Y(Z, t ) ] U 0 ’ ’ (Y) 

+ Z 2 u 23 (Y) - XA(X,Z,t)[V 1 ”(Y) + w 1 ’(Y)]} + .... 

(2.14) 

V = Vj(Y) - A x U 0 (Y) 

+ 5 (XV 2 (Y) - AA X U 0 ’ (Y) 

+ (XA) x [Vj ’ (Y) + wj (Y) ] 

- (ZA) z wj (Y) ) } + . . . , (2.15) 

w = 5[Zwj (Y) ] 

+ 5 2 (XZw 2 (Y) + ZA ( X , Z , t ) W] ' (Y) ) + (2.16) 

P = P 0 + 5 Pj X + i 5 2 P 2 X 2 

+ 1 5 2 P 3 z 2 + 5 5 ^ 2 P(X,Z, t) + ... (2.17) 

The central feature here is the inclusion of the arbitrary function 
y(Z,t), which replaces the arbitrary function of Z (only) found in 
I. In order that a proper match is made with the non* interact ing zone, 
we must have that 

A(X,Z, t ) - p(X 2 X 2 +Z 2 )i 

as X 2 + Z 2 -» ». (2.18) 

The solution (2.14) docs not satisfy the no-slip condition on 
Y = 0, and so it is necessary to include a further (viscous) layer 
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(2.19) 


(•'lower deck’’). This is described by the layer wherein 
Y = Y 5- 1 / 4 = 0(1) 


(see I , for example) . 

The solution then develops in a manner similar to (2.14) - (2.17), 
except (i) Y is replaced by 5 1 / 4 Y, and (ii) there is a correction 
to the u expansion of 5 2 U(X,Y,Z,t), to the V expansion of 
5^/ 4 7(X,7,Z,t), and to the w expansion of 5 2 W(X.Y,Z,t). 

The governing equations for these correction terms are then 


^ + + = 0. 

ax a? az 


( 2 . 20 ) 


s =2 V 3T + 1 n V 2 * a 2 T 0 




dP t a 2 0 
ax a? 2 ’ 

. n 5w a 2 w ap 

i a 2 Y z — = — - — . 

ax aY 2 az 


( 2 . 21 ) 


( 2 . 22 ) 


where a 2 = U o ’’(0). (2.23) 

The boundary conditions to be applied to this system are 

U, W -♦ 0 as Y -» «, (2.24) 

W(Y=0) = V(Y=0) = 0, (2.25) 


and 

0(Y=O) = - I 2 (A 2 +y-p 2 (X 2 X 2 +Z 2 )} . (2.26) 

This system is similar to that found in I, but with the inclusion 
of the time derivative term in (2.21). We now seek to determine the 
condition to be satisfied by A(X,Z,t) if the above system ((2.20) - 
(2.26)) has a solution. 

Differentiating (2.21) with respect to X, (2.22) with respect to Z 
and adding yields 
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, v 2 r a 2 u . a 2 w i . a « av 
* a 2 Y [ ax* + axazj + 32 ? ax 


(2.27) 


a 2 p , a 2 p . a 3 u , a 2 w „ v c a 2 A 

= - + + + • aj I o. . 

ax 2 az 2 ayax ayaz axat 

Invoking continuity, namely (2.20), and differentiating (2.27) with 
respect to Y yields 


d 4 V , v2 d 3 V . _ av 

av 4 a? 2 ax ax 

o a 2 A 

- • s 32 axar 

This equation must be solved subject to 


(2.28) 


V -♦ 0 as Y -» », (2.29) 

whi 1 st on Y = 0 

V = 0, 


— = i a 2 [2AA x -2p 2 X 2 X] , 

a ? 

ill = . [ ill + — 1 . (2.30) 

a ?3 1 ax 2 az 2 J 

Using the analyses of Stewartson (1970), Smith (1982) and I, we 
have that a solution to this system is only possible if the following 
condition is satisfied 

i a 2 [A 2 +y(Z, t) - p 2 (X 2 X 2 +Z 2 )] 


= - (-i)l 


a 2 i 2 3 / 2 (i!) 


dP * f 52P HP 

a^ + J az* 


(X-$)i 




X 

23/ 4 a 9 iS f 3A d^ 

(- 3 / 4 )! J ar ix^tt- 

- oo 


(2.31) 
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The system is closed by recourse to the upper deck, where 
Y = 0(5~3/2) aru j this yields (Smith et al 1977, Duck and Burggraf 1986) 


P(X,Z,t) 



9 2 A/9£ 2 d£dC 

f(X-^) 2 +(Z-C)]i 


(2.32) 


The combined system of (2.31) together with (2.32) represents a closed 
problem, (for prescribed y(Z . t ) ) . We discuss the numerical solution to 
thi s problem next . 
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3- Numeric al Method 


Following I, we choose to use a (double) Fourier Transform method 
to solve the system (2.31) - (2.32). If we define the quantity (for 
example) 

oo oo 

A**(k.ft.t)=J | A(X.Z,t)e-' kX - i,lZ dXdZ. (3.1) 


then in (k.G) space, (2.31) : (2.32) may be easily combined to yield 
(after normalisation to remove positive constants) 

[A 2 +y(Z.O - (X 2 X 2 +Z 2 )]** 

= - ( ik)i (k 2 +H 2 )U** - ( i k) -3/4 S A t **. (3.2) 

The numerical technique used in I was based on assuming the right - 
hand-side of (3.2) (without, of course the time derivative term) was 
known at each iteration level, and then an appropriate method was used 
to solve for the l e f t -hand- s i de terms. Here, we anticipate using a 
time marching scheme to handle t ime -var i a t ions ; however such a scheme 
is likely to require the time-derivative term to be included on the 
left-hand-side of (3.2), otherwise numerical instabilities are very likely 
to occur. Consequently, we choose to transfer the lime-derivative term 
on to the left-hand-side of (3.2), and, in effect to express this term in 


physical (X,Z) space, yielding 


[ A 2 +y(Z, t) - (X 2 X 2 + Z 2 ) 


S 

+ mw 


X 


I 


- oo 


A,(£,Z,t )d5 


] 


= -(ik)i (k 2 +ft 2 )i A** . (3.3) 

Further, just as in I, since A(X,Z,t) is unbounded as 
X 2 + Z 2 -* «, this function is inappropriate for numerical calculation 
(in particular for a Fourier Transform method) in its present form. 

We follow I, and write A(X,Z,l) in the following form 
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(3.4) 


A (X,Z,t) = p 


Y(Z. t) + a n 2 

2p 


+ A r 


+ B, 


where p = (A.2x2+z2+a 0 ^)i , (3.5) 

<2A 0 p)“ = -<ik)i (k2«2> { p = H'(Z . O-°n iL }" (3.6) 

and (Xq is a numerical constant which we have (art i fical ly) introduced 
into the problem, to ensure boundedness as x2 + Z 2 -* 0 (indeed, in 
the limit of zero truncation error, the solution for A(X,Z,t) is 
independent of a 0 . 

We treat B(X,Z,l) in (3.4) as our unknown, described by the equation 

r [y(Z, t )+a Q 2 ]£ + 2Bp 
1 4p2 

- fY(Z.t) + a 0 2] [ ^ | J + A 0 2 

X ** 

+ 2 A 0 B + B2 + } 

= -(ik)l(k’ -f-CL ’ ) 4 (A 0 +B)**. (3.7) 


A ( fu 1 1 y - i mp licit) Crank-Nicol son scheme was applied to (3.7) with 
the initial conditions (at t = 0) prescribed by the corresponding steady 
solution for y(Z, t = 0) , using the computer code developed for I. 

We restricted attention to forcing functions of the form 

Y(Z,t) = Yl(Z) Y2(D. (3- g ) 

which enabled the various "invariant" functions (e.g. A 0 ) to be 
evaluated just once, prior to the unsteady computation commencing. 

One additional feature here, not present in I, is the inclusion of 
the time derivative term. At each timestep (at the ith X location) we 
make the following approximation 
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Xj 

r B.(£,Z,t-At/2)d5 

J -%tr — 

- oo 

Xj -AX 

~ B [ ,Z , t -— w) 

L 

(Xj-^M 

- OO 

+ (AX)3/ 4 { ^ B t (X t ,Z, t - j-) 

- ijyj B( (Xj+AX.Z, t - ) 

+ ^ B t (Xi-AX,Z,t-jl) }, (3.9) 

where AX is the grid size in the X direction (and AZ the grid size 
in the Z direction), and At the time step. The Crank-Ni col son 
diffencing scheme made the following approximation 
B t (Xj ,Z, t - j~) 

B(X; ,Z,t)-B(Xi ,Z, t -At ) (310) 

_ 

whilst the integral term in (3.19) was evaluated using a standard 
trapezoidal scheme (coupled with (3.10)). 

The term involving B on the right-hand-side of (3.7) was evaluated 
using the fast Fourier Transform technique of Cooley and Tukcy (1965), 
in which the function is actually evaluated in spectral ( k , 0. ) space, 
and then transformed back into physical (X,Z) space; this same 
technique was employed in I, and results in considerable savings on 
computational times. 

At each time and iteration level, at each Z location, the (algebraic) 
approximation to (3.7) can be written symbolically as 

I B 2 + E B = R, (3.11) 

where I is the unit matrix and B = B(Xj(i=l,N), Z.t). 


We then write 
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( 3 . 12 ) 


B = B( n_1 ) + 5B 

where B( n ) denotes the value at the n th iteration level; here we 
discard 0(|5B| 2 ) terms, to yield the system 
(2 l B(") + E )5B 

= R - I B(n) 2 - E B< n > (3.13) 

Thus to obtain a new estimate for the B(°) a linear algebraic system was 
solved (in I, just a scalar system was solved; the off-diagonal terms in 
(3.13) arise directly from the time derivative term.) This scheme is 
loosely based on Newton’s method. The solution was then determined at all 
Z stations, and the overall procedure then repeated until convergence 
was attained, where-upon the calculation proceeded to the following 
t imestep. 
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4. Numerical Solutions 


In this section we present numerical solutions for two choices of 
y(Z,t) (and, as it turns out these exhibit quite different features). 


The first example 

taken 

was with \ = S = 1 and 


Y(Z,t) = 

-Z 2 
e L 

(-2+4tanh 2 t) , for t > 0, 

(4.1) 

and hence we have 




Yl (Z) = e 

•z 2 , 


(4.2) 

Y2O) = - 

2 + 4 

tanh 2 t , t > 0. 

(4.3) 


In this example (and in the following example) conditions at l = 0 
were obtained by solving the steady problem (considered in I) with 
Y(Z,0). Further, since Y 2 O) varies (monotonical ly) from -2 at 
t = 0 to +2 as t -+ «*, then according to I, the corresponding steady 
problem with (4.2) lies entirely within the envelope of steady solutions 
found in I. 

Figure 1 shows the temporal variation of A(X,Z,t) at selected X 
and Z stations, and these confirm an evolution from the initial steady 
state Yl ( t=0)=-2 to the final y 1 C t » I ) = + 2 state (and Fig.l 
shows results agreeing with the corresponding steady "upper branch" 
solution in this limit). Interestingly, a common trend observed at the 
locations considered in Fig.l is that, although Y(Z.t) reduces 
monotonical ly with time (as would quasi-steady results for A(X,Z), based 
on the results of I), there is seen to be a very slight "tough in 
the results, soon after the computation commences, followed by a very 
pronounced peak soon after t = 1, followed by a monotonic decrease towards 
the t » 1 values. This particular computation was carried out using 
a grid which extended from X (and Z) - 10.32 to X (and Z) — 9.68 
with AX = AZ = 0.64, At = 0.005 and a 0 = I ; we shall refer to this 
as grid I . 

The second (and what turns out to be the more challenging and 
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interesting) example taken was X = S = 1 with 

Y(Z,t) = e‘ z2 (-2+8tanh 2 t) , t >0, (4.4) 

and hence y ^ (Z) is given by (4.2), whilst 

Y 2 ( t ) = - 2+8tanh 2 t, t > 0. (4.5) 

Thus we see 72(0 extends from -2 at t = 0 (as in the previous 
example) to +6 as t However this latter value lies outside 

of the envelope of steady solutions found in 1 corresponding to our 
choice of Yl (Z) . 

The distribution of A(X,Z) at fixed X or Z stations and fixed 
t arc shown in Figs. 2 - 5. These particular computations were performed 
on a grid that extended from X (and Z) = -10.15 to X (and Z) = 9.84, 
with AX = AZ = 0.32, At = 0.005 and a 0 = 1; we refer to this as grid 
II. 

It is quite clear that a singularity is forming, close to X = Z = 0, 
at a finite time, with A(X,Z,t) becoming progressively larger and more 
negative at this location. Figure 6 shows the temporal variation of 
A(X=Z=0) on a variety of grids, to enable a (partial) assessment to be made 
of the accuracy of the scheme in this critical region. The solidus denotes 
the grid II results. The broken line - - - represents the grid I results 
(a computation was also carried out on this grid, but with a Q = 3; 
these proved to be indi st igui shable on the scale used in Fig. 6). 

The dot-dashed line denotes grid III results with X 

and Z extending from (approxima t ly) -5.16 to + 4.84 with 

AX = AZ a 0.32, At = 0.005 and a 0 = 1 . Grid IV results are shown as a 

dotted line a grid which is identical to grid I. except 

At = 0.05. All these results confirm the likely appearance of a finite-time 
singularity close to X = Z = 0; the accuracy of the computations is also 
con f i rmed . 

Figure 7 shows the (spectral) distribution of Re j B**(k ,d=0, t ) j 


16 



with k at selected times, (computed on grid II) It appears that this 
distribution (and indeed other spectral distributions, not shown here) 
undergoes a rapid change in behaviour as the apparent breakdown is 
approached. In particular the decay of the solution as iki » 1 is seen 
to diminish (and this in turn mirrors the finite-lime, f ini te- local ion 
breakdown in physical X,Z space). 

In the following section. we go on to analyse the structure of this 
likely si ngu lari t y . 
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5 . The breakdown 


We now seek to determine the nature of the breakdown suggested in a 
number of the numerical results described in the previous section. 

In particular, we seek an asymptotic solution which predicts 
I A i -♦ « at a finite value of X (X 0 say), at a finite value of Z 
(Z D say); (notice that the indications from the results for the 
configuration in the previous section that resulted in a breakdown are 
that Z 0 = 0), at a finite line (t s say). We define the (small) 
t imescale 

x = t s - t; (5.1) 

then inspection of (3.7) suggests that if the (i) non-linear, (ii) 
time derivative and (iii) the right-hand-side terms all balance, then 
the following (‘‘hatted’’) variables are the most appropriate 
X - X 0 = x 4 / 9 S' 4 / 9 X 
Z - Z 0 = T 4 / 9 S' 4 / 9 z 

A = I" 2 / 3 S 2 /3 A (X.Z) + o(x-2/3). (5.2) 

Substitution of these forms into (3.3), taking the largest terms 
in x, and writing the solution entirely in physical (X.Z) space yields 
the following system 


A 2 


: A + 


1 *i + 1 A c 4 


9n») 

4 


(X - 5)* 



oo 


- oo 


oo 

.. ^ i i + . A c i i 1 K(x> dwe, 

[(X-5)2+(Z-02]i 

- 00 


(5.3) 


where 


*• ( 


1 + 


X-^ 

[(X-^)2+(Z-0 2 ]i 




(5.4) 
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and K(x) denotes the complete elliptic integral of the first kind, 
argument %. 

In order to match correctly on to the X = 0(1), Z = 0(1) problem, 
we requi re that 

A - a(X/Z) (X!X 2 +Z 2 )-3/4, (5.5) 

A ^ A A 

as X 2 + Z 2 -» », where Xj is a constant and a(X/Z) is a function of 
X/Z. Both quantities must be determined from the terminal behaviour of 
(3.2). (The scalings (5.2) arc entirely consistent with the two- 
dimensional work of Smith 1982). 

Indeed, it is possible to obtain a result identical to (5.5) using 
a linearized form of (3.2). This may be justified by supposing that the 
parameter S » 1, i.e. we arc concerned with rapid temporal variations in 
external forcing. We also suppose that the typical X and Z lcngthscales 
are both 0(S' 4 / 9 ). The resulting system is then (having scaled out S) 

A** + (k 2 +ft 2 )i ( i k)5/4 A ** = 0 . (5.6) 

The solution to this is 

A** = A 0 **(k,ft) e(k 2 +fi 2 ) i (*k) 5/4 0-t s ). (5.7) 

From this we see that as t -* t s , then the solution approaches a finite- 

4 

time breakdown, with the appropriate X and Z scales being 0((t s -t) 9 ), 
in both cases, in complete accord with (5.2). 

Let us now return to consider the system (5.3); this, in general, 
requires a fully numerical (and quite substantial) task. However, we 
can make some further analytic progress (that suggests the existence of 
regular solutions to (5.3) for all X and Z) , by using the idea, 
put forward by Smith (1982), of assuming |A| « 1 for all X and Z (this 
amounts to assuming i a (X/Z) i « 1). With this restriction the ‘A 2 ' 

term on the 1 c f t -hand - s i de of (5.3) may be discarded. If the resulting 
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system is subject to a double Fourier transform, then (after 


normalising to remove positive constants, and then using tilde variables) 


j A** + k + ft 


9k 


0a 


= (ik)5/4 (k2+B;2)i a** . 


(5.8) 


This may be solved using a technique based on the method of 

~2 2~2 i ~ - 

characteristics; defining variables £ \ = (k +\jft )* , £2 = 

(5.8) may then be written in the form 


$1 


0A** 

^1 


{- 


1 + 1 L } A* 

(1+ — )5/8 ^2+Xj2 

^2 


(5.9) 


This can be integrated routinely to give the following general solution 

A**(k,B) = F**( § )(k 2 +X 1 2 a 2 )" i exp {-^(ik) 5 / 4 (k 2 +ft 2 ) 1 }. 

* (5.10) 


Here F**(k/ft) is an arbitrary function of (k/ft). The solution (5.10) is 
entirely consistant with (5.5) (seen by allowing k 2 + ft 2 -» 0). The 
importance of (5.10) is that it illustrates that a regular solution 
exists to (the linearised form of) (5.3) for all X and Z, and as 
such strongly suggests that the terminal description of the breakdown 
is given by the solution of (5.3). 
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6. Genera l Discussion 

In this paper (i) we have presented the structure of the flow in 
the vicinity of a three-dimensional marginal separation point (situated 
along a line of symmetry), including the effects of unsteadiness at 
this point; (ii) the solution of the non-linear, unsteady integro- 
diffcrential system has been carried out in two distinct cases, the 
first where the flow evolves regularly from one steady state to another 
steady state; and the second where a three-dimensional breakdown is seen to 
be observed; (iii) the structure of this breakdown has also been described. 

Although, from the point of view of this paper, the particular 
details of the mechanism by which unsteadiness is injected into the 
problem are not important, we could, perhaps envisage this may be 
achieved (for example) by a pitching/buffeting of an aerofoil/ or a 
aerodynamic body. 

Qualitative comparison between the numerical results of Section 4, 
with the analysis of Section 5 seems quite supportive, although the 
computations do become very much more difficult as the breakdown is 
approached (see the remarks of Smith 1982). Certainly the numerical 
results point to the emergence of a small region (close to X = Z =0) 
inside of which A(X,Z,t) increases rapidly, both spatially and temporally. 

The implications of this breakdown are that as Y( Z '0 increases 
above its critical value although the flow initially responds slowly 

I 

(on a timescale 0(R*M), the flow then develops sudden 1 y/explos i vc 1 y as 
the breakdown time is approached. 

A further matter concerns the ulimate behaviour of the flow, 
leading from the breakdown detailed in Section 5. It is surmised by 
Smith (1982) that (in the two-dimensional context) the breakdown 
leads ultimately to the flow being governed by a form of the unsteady 
boundary- layer equations. Here, in the three-dimensional situation, such 
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a scenario seems not unreasonable. This aspect is currently under 
invest igat ion. 
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